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Abstract  _ _ 

Calculations  are  undertaken  for  the  electronic  partition  functions  of  plasmas  in  a 
temperature  and  pressure  range  relevant  to  electrothermal-chemical  gun  applications. 
The  intent  is  to  assess  the  limitations  of  results  obtained  previously,  as  well  as  to  provide 
an  improvement  for  existing  models  for  the  capillary  discharge.  Various  thermophysical 
properties  that  depend  on  the  electronic  partition  functions  are  calculated  and  compared 
with  results  obtained  in  preceding  work.  The  general  conclusion  reached  is  that  accurate 
predictions  of  the  concentrations  of  minority  species  in  the  plasma  require  highly 
accurate  calculations  of  the  partition  functions.  Less  accurate  approximations,  however, 
predict  the  concentration  of  majority  species,  as  well  as  most  of  the  thermophysical 
properties  of  the  plasma,  reasonably  accurately. 
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1.  Introduction 


The  purpose  of  this  note  is  to  undertake  calculations  for  the  electronic  partition  functions  of 
plasmas  composed  of  the  constituents  of  polyethylene.  Since  the  calculations  are  more  accurate 
than  those  of  our  previous  work,  the  analysis  is  a  means  of  evaluating  the  limitations  of  existing 
models  and/or  extending  their  applicability  to  a  wider  range  of  temperatures  and  pressures.  Only 
polyethylene  is  considered  because  of  its  wide  use  as  a  plasma-producing  liner  in  capillaries 
designed  for  electrothermal-chemical  (ETC)  launchers.  However,  it  is  expected  that  the  general 
conclusions  of  this  note  should  apply  qualitatively  to  other  substances  as  well. 

In  addition  to  the  general  desire  to  provide  calculations  that  are  as  accurate  and  internally 
consistent  as  possible,  undertaking  this  calculation  has  also  been  motivated  by  some  recent  work 
performed  by  McQuaid  and  Nusca  (1999).  In  that  work,  the  authors  noted  that  while  the 
approximations  employed  in  our  earlier  work  seemed  adequate  for  the  cases  of  interest  then,  they 
might  prove  inadequate  under  certain  other  conditions  (the  details  of  which  will  be  discussed 
subsequently).  The  intention  of  the  present  work  is  two-fold.  The  first  goal  is  to  indicate  what 
specific  plasma  properties  are  not  adequately  predicted  with  the  approximations  in  the  existing 
models  and  under  what  conditions  that  failure  occurs.  The  second  goal  is  to  perform  more 
accurate  calculations  and  include  the  methodology  for  general  use  in  the  existing  models.  We 
will  refer  extensively  to  our  previous  work  (Powell  and  Zielinski  1992,  1993)  in  order  to  avoid 
needless  repetition.  The  first  of  these  references  is  referred  to  as  PZ1 . 

Various  electrical  and  thermodynamic  properties  of  a  plasma  depend  on  the  electronic 
partition  functions  of  its  constituents.  This  function  for  any  particular  monatomic  constituent  a 
can  be  written  as  (see  equation  [24],  PZ1) 

Zjt.=Zgi.iexp(-Uj„,/kT),  (1) 
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where  Ujaj  is  the  energy  level  of  the  ith  electronic  state  of  the  jth  ion,  and  gjaj  is  the  appropriate 
degeneracy  factor  for  this  level.  The  actual  calculation  in  equation  (1)  is  problematic  for  several 
reasons.  First,  for  an  isolated  atom  or  ion,  the  sum  diverges.  This  divergence  arises  because 
there  are  an  infinite  number  of  discrete  states,  and  the  levels  Ujaj  approach  a  constant  value, 
namely,  the  ionization  energy,  as  i  becomes  large.  This  problem  is  somewhat  illusory  since  the 
radii  of  the  electronic  orbits  increase  with  increasing  energy,  eventually  becoming  comparable  to 
the  interatomic  spacing.  At  this  point,  the  atom  is  essentially  ionized.  The  effect  is  generally 
modeled  by  reducing  the  ionization  potential  of  the  isolated  atom  or  ion  by  an  amount  that 
depends  on  the  density  of  the  mixture,  and  by  including  in  the  summation  only  the  terms  whose 
energy  is  less  than  the  reduced  ionization  energy.  The  sum  therefore  remains  finite.  However, 
general  agreement  does  not  exist  on  the  amount  by  which  the  ionization  potential  should  be 
reduced. 

The  second  problem  arises  because  the  energy  levels  in  equation  (1)  are  inadequately  known 
either  experimentally  or  theoretically.  Most  investigators  rely  on  the  data  contained  in  Moore’s 
famous  tables  (Moore  1971),  but  many  of  the  levels,  particularly  those  that  correspond  to  high 
orbital  angular  momentum  quantum  numbers,  have  never  been  observed  experimentally  and  are 
missing  from  the  tables.  Thus,  one  is  forced  to  either  ignore  the  missing  levels  (Kovitya  1985) 
or  estimate  their  values  from  rather  crude  theoretical  considerations. 

Simpler  approximations  are  also  common  in  calculations  involving  the  partition  functions. 
In  some  cases  (Loeb  et  al.  1987;  Mohanti  1990),  the  ratio  of  various  functions  that  occurs  in 
certain  equations  is  replaced  by  unity  as  a  crude  approximation,  and  the  sum  is  not  even 
evaluated.  Other  authors  employ  only  the  first  few  terms  (Cambel  1963;  Zeldovich  and  Raizer 
1966),  or  even  just  the  leading  term  (i.e.,  the  degeneracy  of  the  ground  state)  (Zaghloul  et  al. 
2000;  Ngo  et  al.  1998)  in  the  sum  in  equation  (1).  Presumably,  these  approximations  that 
employ  only  a  few  terms  are  made  with  the  belief  that  the  higher  energy  states  are  not 
appreciably  populated.  While  the  approximations  frequently  seem  to  work,  they  are  difficult  to 
justify  in  a  rigorous  sense  because  the  degeneracy  factors  increase  with  increasing  principal 
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quantum  numbers.  Thus,  although  the  higher  energy  states  are  not  highly  populated,  they  may 
have  a  significant  effect  because  of  their  greater  number. 

In  our  earlier  calculations  of  plasma  properties  for  both  electromagnetic  (Powell  and  Batteh 
1981)  and  electrothermal  (Powell  and  Zielinski  1992,  1993)  launchers,  only  the  first  few  terms 
were  generally  included  in  the  partition-function  calculations.  For  the  polyethylene  plasma,  only 
those  terms  from  Moore’s  tables  that  had  energy  less  than  about  9  ev  were  included.  Some 
studies  were  undertaken  in  which  the  number  of  terms  was  varied  and  only  minimal  changes  in 
the  plasma  properties  were  observed,  but  these  studies  were  not  extensive.  Furthermore, 
McQuaid  and  Nusca  (1999)  have  pointed  out  that  under  some  conditions,  discrepancies  exist 
between  the  magnitude  of  partition  functions  calculated  with  only  those  terms  retained  and  the 
functions  computed  by  Drawin  and  Felenbok  (1965).  The  discrepancy  occurs  at  high 
temperatures  and  at  pressures  generally  lower  than  those  relevant  to  our  calculations.  In  the 
interest  of  being  as  complete  and  thorough  as  possible,  we  sought  to  obtain  better 
approximations  for  the  functions  and  to  determine  their  impact  on  various  plasma  properties. 

2.  Calculations 

As  in  PZ1,  the  plasma  is  assumed  to  be  composed  completely  of  the  monatomic  constituents: 
H,  H+,  C,  C+,  C++,  and  e~.  For  polyethylene,  there  are  twice  as  many  hydrogen  heavy  particles 
(ions  and  neutrals)  as  carbon  heavy  particles.  As  before,  we  also  assume  that  the  ionization 
potential  for  the  jth  ion,  Ija,  is  reduced  because  of  the  finite  density  by  an  amount  AIja  (Ebeling 
and  Sandig  1973)  given  by 


AT,= - - = - T. 

J  4tts0(Xd  +  A  /8) 


(2) 


In  this  equation,  e  represents  the  electron  charge  and  so  the  permittivity  of  free  space.  The 
parameters  Ad  and  A  are  the  Debye  length  and  de  Broglie  wavelength,  respectively.  These 
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quantities  depend  on  the  temperature  and  pressure  of  the  plasma;  therefore,  they  must  be 
calculated  self-consistently  with  the  other  plasma  properties.  Explicit  formulas  for  their 
calculation  are  given  in  equations  (20)  and  (21)  in  PZ1.  In  contrast  to  previous  calculations  in 
which  the  sum  in  equation  (1)  was  restricted,  included  in  this  calculation  are  all  terms  whose 
energy  is  less  than  the  reduced  ionization  potential  (i.e.,  Ija-  AIja). 

For  hydrogen,  all  of  the  relevant  states  needed  have  been  calculated,  and  they  are  in  Moore’s 
(1971)  tables.  Thus,  computation  of  the  function  Z0H  poses  no  problem.  For  carbon,  on  the 

other  hand,  Moore  lists  levels  that  have  been  observed  experimentally,  and  the  tables  are 
incomplete.  In  particular,  some  levels  corresponding  to  high  energy  states  had  not  been  observed 
at  the  time  the  tables  were  compiled,  and  those  levels  are  missing.  To  rectify  that  situation, 
McQuaid  and  Nusca  (1999)  attempted  to  fill  in  the  missing  levels  both  with  experimental  data 
from  Stiganov  and  Sventitskii  (1968),  as  well  as  with  their  own  estimates.  Their  table,  compiled 
for  neutral  carbon,  is  employed  in  the  current  calculation  to  compute  the  function  Zoc .  For  ZIC 
and  Z2C  (i.e.,  the  functions  that  correspond  to  first-  and  second-ionized  carbon),  only  the  terms 

that  have  been  observed  and  are  included  in  Moore’s  tables  have  been  accounted  for.  McQuaid 
and  Nusca  (1999)  employed  a  similar  approximation,  actually  including  those  levels 
corresponding  to  an  energy  less  than  about  1 1  ev.  They  have  reasonably  argued  that  the  higher 
energy  states  should  contribute  negligibly  in  the  temperature  range  of  interest.  However,  in  the 
present  work  the  higher  levels  were  included  anyway,  but  the  sums  were  truncated  at  I!C  —  AI1C 

and  I2C  -  AI2C .  All  of  the  levels  in  the  calculations  for  the  various  functions  are  indicated  in  the 
Appendix. 

Calculations  of  all  plasma  properties  have  been  performed  for  temperatures  ranging  from  1  to 
3.5  ev  and  pressures  ranging  from  1  to  400  MPa.  Quantities  calculated  for  each  temperature  and 
pressure  include:  (1)  the  internal  energy  E,  (2)  the  density  p,  and  (3)  the  ionization  factors  xtH, 
xic,  and  x2c-  These  factors  represent  the  ratio  of  first  or  second  ions  for  the  particular  species  to 
the  total  number  of  heavy  particles  for  that  species.  All  of  these  quantities  can  be  calculated 
using  equations  (15)-(23),  as  well  as  equation  (25),  of  PZ1.  It  is  necessary,  however,  to  solve 
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for  the  various  properties  iteratively,  since  the  nonlinear,  algebraic  equations  that  determine  their 
values  are  coupled.  Only  results  are  presented,  and  the  reader  is  referred  to  PZ1,  in  which  a 
similar,  but  more  complicated,  calculation  is  undertaken  for  details. 

3.  Results  and  Discussion 

In  Figures  1  and  2,  E  and  p  are  plotted  as  a  function  of  the  temperature  and  pressure,  and  in 
Figures  3-5  are  the  ionization  factors  xih,  xjc,  and  X20  These  particular  variables  were  selected 
for  plotting  since  most  other  quantities  of  interest,  such  as  the  electrical  conductivity  or  electron 
density,  follow  in  a  very  simple  way.  These  quantities  are  plotted  as  a  function  of  temperature  T 
at  four  different  pressures  P.  The  solid  curves  correspond  to  results  obtained  when  only  levels 
less  than  9  ev  were  included  in  the  partition  functions,  as  was  done  in  the  original  calculations; 
the  dotted  curves  correspond  to  results  obtained  when  all  levels  less  than  the  reduced  ionization 
potential  were  included  in  the  manner  described  previously. 

The  results  of  the  calculations  exhibit  fairly  obvious  behavior.  The  internal  energy  increases 
with  increasing  temperature  at  a  fixed  pressure,  and  it  decreases  with  increasing  pressure  at  a 
fixed  temperature.  The  latter  result  follows  because  as  the  pressure  increases,  the  interatomic 
spacing  becomes  smaller,  and  collisions  that  tend  to  deionize  the  gas  become  more  prevalent. 
Consequently,  at  low  pressures  a  substantial  part  of  the  energy  of  the  plasma  is  associated  with 
the  energy  required  to  produce  ionization.  In  a  similar  manner,  the  ionization  factors  xih  exhibit 
the  same  type  of  behavior.  For  carbon,  similar  behavior  is  observed  at  lower  temperatures,  but 
eventually  second  ionization  becomes  important,  and  the  relative  number  of  first  ions  is  reduced. 
At  the  high  temperatures  and  very  low  pressures  considered,  particularly  at  P=1  MPa  with  T 
greater  than  about  2.5  ev,  third  ionization  is  probably  important,  but  it  is  not  accounted  for  in  the 
model.  However,  such  low  pressures  are  generally  outside  the  range  of  interest  for 
electrothermal  or  electromagnetic  gun  applications.  The  mass  density  decreases  with  increasing 
temperature  at  constant  pressure  to  satisfy  the  equation  of  state. 
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Figure  1.  Internal  Energy  as  a  Function  of  Temperature  and  Pressure.  Solid  curves 
represent  results  when  only  levels  less  than  9  ev  are  retained  in  the  calculations; 
dotted  curves  represent  results  when  all  terms  less  than  the  reduced  ionization 
potential  are  retained.  Pressures  are  denoted  by:  black,  P  =  1  MPa; 
red,  P  =  10  MPa;  blue,  P  =  100  MPa;  and  green,  P  =  400  MPa. 


Figure  2.  Density  as  a  Function  of  Temperature  and  Pressure.  Solid  curves  represent 
results  when  only  levels  less  than  9  ev  are  retained  in  the  calculations;  dotted 
curves  represent  results  when  all  terms  less  than  the  reduced  ionization 
potential  are  retained.  Pressures  are  denoted  by:  black,  P  =  1  MPa;  red,  P  = 
10  MPa;  blue,  P  =  100  MPa;  and  green,  P  =  400  MPa. 
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Figure  3.  First  Ion  Concentration  of  Hydrogen  as  a  Function  of  Temperature  and 
Pressure.  Solid  curves  represent  results  when  only  levels  less  than  9  ev  are 
retained  in  the  calculations;  dotted  curves  represent  results  when  all  terms  less 
than  the  reduced  ionization  potential  are  retained.  Pressures  are  denoted  by: 
black,  P  =  1  MPa;  red,  P  =  10  MPa;  blue,  P  =  100  MPa;  and  green, 
P  =  400  MPa. 


Figure  4.  First  Ion  Concentration  of  Carbon  as  a  Function  of  Temperature  and  Pressure. 

Solid  curves  represent  results  when  only  levels  less  than  9  ev  are  retained  in  the 
calculations;  dotted  curves  represent  results  when  all  terms  less  than  the  reduced 
ionization  potential  are  retained.  Pressures  are  denoted  by: 
black,  P  =  1  MPa;  red,  P  =  10  MPa;  blue,  P  =  100  MPa;  and  green,  P  =  400  MPa. 
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Figure  5.  Second  Ion  Concentration  of  Carbon  as  a  Function  of  Temperature  and 
Pressure.  Solid  curves  represent  results  when  only  levels  less  than  9  ev  are 
retained  in  the  calculations;  solid  curves  represent  results  when  all  terms  less 
than  the  reduced  ionization  potential  are  retained.  Pressures  are  denoted  by: 
black,  P  =  1  MPa;  red,  P  =  10  MPa;  blue,  P  =  100  MPa;  and  green, 
P  =  400  MPa. 

Of  the  greatest  interest  in  the  results  is  the  very  weak  dependence  on  the  number  of  terms 
included  in  the  partition  functions.  In  most  cases,  differences  are  imperceptible.  For  example, 
consider  the  results  for  xm  at  P  =  100  MPa  and  T  =  3  ev,  a  point  where  one  of  the  greater 
deviations  is  observed.  Superscript  s  (for  solid)  denotes  results  for  the  original  calculation,  and 
superscript  d  (for  dotted)  denotes  results  for  the  case  in  which  all  terms  are  included  in  the 
partition  functions.  At  the  temperature  and  pressure  in  question  we  find  x'H  =  0.7555  and 
xfH  =  0.7076,  a  difference  of  only  about  7%.  On  the  other  hand,  the  partition  functions  at  this 
point  (not  plotted),  given  by  ZS0H  =2.000  and  Z„H  =  2.604,  differ  by  23%.  The  reason  that 

significant  differences  in  the  partition  functions  for  the  two  calculations  lead  to  only  small 
differences  in  the  relative  number  of  ions  can  be  described  as  follows.  At  temperatures 
sufficiently  high  that  there  is  appreciable  population  of  the  higher-energy  states  (i.e.,  the  terms 
omitted  in  the  calculation  of  Zsm  should  not  really  have  been  neglected),  most  of  the  hydrogen 
atoms  (71%)  have  already  been  ionized.  Thus,  if  the  temperature  is  increased  still  more  so  that 


8 


ZgH  is  even  larger  relative  to  Zsm ,  minimal  changes  will  be  observed  in  xfH  since  its  value  can 
never  exceed  unity.  To  make  the  argument  more  concrete,  consider  the  values  obtained  at  the 
same  pressure  and  at  T  =  3.5  ev,  where  x*H  =0.8720,  ZS0H  =2.000,  xfH  =0.8222,  and 
Zdm  =  2.990.  For  this  case,  the  percentage  differences  in  the  partition  functions  are  33%, 

whereas  the  differences  in  the  ionization  factors  are  only  6%.  Similar  conclusions  apparently 
hold  (perhaps  even  more  so)  for  carbon,  except  that  xi  approaches  zero  rather  than  unity  at  high 
temperatures  when  the  second  ionization  begins  to  dominate. 

It  is  important  to  note  that  while  these  conclusions  hold  for  the  majority  species,  they  do  not 
hold  for  the  minority  species.  Thus,  the  concentrations  of  neutral  hydrogen  predicted  at  T  =  3  ev 
are  XqH  =0.2445  and  x^H  =0.2924;  at  3.5  ev  the  corresponding  values  are  XqH  =0.1280  and 

XqH  =  0.1778.  Here,  the  percentage  differences  are  16%  and  28%,  respectively,  and  the  situation 
gets  worse  as  the  temperature  increases. 

Rolader  (1987)  reached  similar  conclusions  in  a  careful  study  of  partition  function 
calculations  for  aluminum  and  copper,  two  substances  that  are  important  in  plasma  armatures  for 
electromagnetic  guns.  In  particular,  he  considered  three  different  methods  of  performing  the 
calculations.  In  the  first,  only  levels  given  in  Moore’s  tables  were  considered;  in  the  second  and 
third,  two  different  approximations  (hydrogenic  and  Rydberg-Ritz)  were  employed  to  fill  in  the 
missing  levels.  In  all  cases,  the  partition  function  series  were  terminated  at  the  reduced 
ionization  potentials.  At  high  temperatures,  where  the  higher-energy  states  were  most  populated, 
and  at  low  pressures,  where  the  reduction  in  ionization  potential  was  small,  the  number  and 
magnitude  of  the  terms  included  in  the  sum  varied  greatly.  The  final  conclusion,  however,  was 
that  in  any  event,  only  the  minority  species  were  significantly  affected.  Fortunately,  most  of  the 
electrical  and  thermodynamic  properties  of  the  plasma  do  not  depend  sensitively  on  the  minority 
species.  It  is  for  this  reason  that  we  can  be  somewhat  cavalier  in  calculating  these  functions. 
Rolader,  however,  pointed  out  that  minority  species  can  be  important  in  calculating  the  opacity 
of  the  plasma,  and  that  spectroscopic  studies  could  focus  on  minority  species  whose  probabilities 


9 


may  be  better  known  than  those  of  the  majority  species.  For  studies  such  as  these,  detailed  and 
careful  calculations  of  the  partition  function  would  clearly  be  required. 
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Table  A-l  contains  the  data  employed  in  calculations  of  the  electronic  partition  functions  for 
H,  C,  C+,  and  C++.  These  data  were  obtained  from  the  sources  indicated  earlier.  In  some  cases, 
several  levels  of  comparable  energy  are  combined  into  a  single  level  denoted  by  i  with 
degeneracy  g.  Maximum  values  included  in  the  tables  are  always  greater  than  the  reduced 
ionization  potential  for  any  of  the  calculations  undertaken. 
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Table  A-l.  Energy-Level  Data  Employed  in  Partition-Function  Calculations  for 
Hydrogen  and  Carbon 
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Table  A-l.  Energy-Level  Data  Employed  in  Partition-Function  Calculations  for 
Hydrogen  and  Carbon  (continued) 
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